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Abstract 

We analyze homogeneous, anisotropic cosmology driven by a self-interacting 
"massive" antisymmetric tensor field which is present in string theories with 
D-branes. Time-dependent magnetic B field existing in the early universe can lead 
to the Bianchi type I universe. Evolutions of such a tensor field are solved exactly 
or numerically in the universe dominated by vacuum energy, radiation, and B field 
itself. The matter-like behavior of the B field (dubbed as "S-matter") ensures that 
the anisotropy disappears at late time and thus becomes unobservable in a reason- 
able cosmological scenario. Such a feature should be contrasted to the cosmology 
of the conventional massless antisymmetric tensor field. 
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1 Introduction 



The antisymmetric tensor field (Kalb-Ramond field) B^ u is an indispensable field in 
effective theories derived from the string theory. It was considered as a massless de- 
gree of freedom in the bulk spacetime and thus it is dual to a pseudo-scalar field h: 
H^ up = d[pB up ] = e fiupa d a h in four dimensions [T]. Once we take into account D-brane 
which easily forms a bound state of Dp-brane-fundamental strings [2], such a B field alone 
becomes gauge non-invariant on the Dp-brane and gauge invariance is restored through 
its coupling to a U(l) vector field living on a Dp-brane, i.e., the gauge-invariant field 
strength is 

+ = B/ju,, where F pv = d p A v - d v A p . (1.1) 

Note that B pv is invariant under the tensor gauge transformation 

SB^ = d^Cu ~ d v Cp , SAp = -( fl /2na' 

including the ordinary U(l) gauge transformation = d^X. Intriguingly enough, the 
effective action of antisymmetric tensor field derived on the D-brane dictates this Abelian 
tensor field of second rank to be massive and self-interacting [HI 121 0] . 

Another indispensable field is gravitational field g pv and, once it also enters into the 
playground in the string scale, it is natural to deal with cosmo logical issues due to the 
existence of the antisymmetric tensor field. In this paper, we consider the cosmology of 
the antisymmetric tensor field in a D-brane world which is supposed to contain our uni- 
verse. The cosmological solutions with the massive antisymmetric tensor field B are quite 
different from those with the massless one. In the latter case, a homogenous and isotropic 
universe can be realized for H 0i j = and H^k = ftijk with some constant / oc d t h [3]. 
Another homogeneous solution leading to an anisotropic Bianchi type I universe takes the 
configuration of B^ = and H 0i j = e^f-Lf. with one non- vanishing constant oc d^h 0. 
Such an anisotropy could be in conflict with observation. Indeed, the anisotropy driven 
by a non- vanishing grows in an expanding universe 0. In relation with anisotropy in 
Bianchi type universes, various issues have been discussed [7]. In the brane world where 
the antisymmetric tensor field becomes massive, time-dependent homogeneous solutions 
are described by "magnetic" degrees of freedom By = Bij{t) leading to H%jk = and 
Hoij = dtBij, which does not allow isotropic cosmology. Assuming the simplest form of 
anisotropic universe, Bianchi type I, we will analyze the cosmology with the massive mag- 
netic field B{j in various situations of inflationary, radiation dominated, or £>-dominated 
era. As we will see, the initial anisotropy is damped away and thus a practically isotropic 
universe is recovered at late time due to the ordinary matter-like behavior of the B field. 

The paper is organized as follows. In Section 2, we describe the four-dimensional 
low-energy effective action of D-brane world from which the mass of the antisymmetric 
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tensor field and its equations of motion coupled to the gravity are derived. In Section 3, 
we find time-dependent homogenous solutions with one magnetic component $12 (t) 7^ 
which leads to the Bianchi type I cosmology. Section 4 concludes the paper. 



2 Antisymmetric Tensor Field on D-brane 

Suppose that our universe is a part of a Dp-brane. Then its bosonic sector comprises bulk 
degrees including the graviton g^, the dilaton $, and the antisymmetric tensor field 
of rank-two in addition to U(l) gauge field living on the Dp-brane. Combining B^ v 
with as in (jl.lj) . the bosonic action of our system in string frame is given by the sum 
of the bulk terms in ten dimensions; 



S 
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2 K io 



d 10 x 



-2$ 



R-2A 10 + 4(d$) 2 --^H 2 



(2.2) 



and the brane terms in p spatial dimensions; 



S Dp = -fM p J d p+ WV-det(<7 + £) 



(2.3) 



\{2-n) 7 a' 4: and fi 2 = (n/ nl )(4:7r 2 a') 3 p 4j. Recall that a' defines the string 
Here the spacetime indices are not explicitly expressed and the field 



1-1/2 



where k\ 
scale; m s = a 

strength H of the antisymmetric tensor field is again H^p = d^B vp \ because of Bianchi 
identity of the U(l) gauge field. 

To consider the cosmology of our universe, let us compactify six extra-dimensions with 
a common radius R c and assume the stabilized dilaton leading to a finite string coupling 
g s = e*. Then, the four-dimensional effective action becomes 



2/?4 



d x\ 



-9 



R-2A-—H fMVp H^ 
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(2.4) 



where B*^ v = \/—gep iVa pB a,3 /2 with 60123 = 1- Here the reduced four-dimensional Planck 
length k 4 is determined by = R®/g 2 Kl . The effective action ()2.4|) up to the quadratic 
terms shows that the antisymmetric tensor field acquires the mass; 



m B 



P-3 
16 



mp 



15-p 
8 



m P , 



(2.5) 



where a s = g 2 /4ir, m c = 1/R C , and mp = 1/k 4 = 2.4 x 10 18 GeV. As we will see, the 
massive B field (named as 5-matter) behaves like a massive scalar in some respects in 
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its cosmological evolution, due to which the intrinsic anisotropy set by directional B field 
gets diminished eventually. 

From the action (|2.4jl . one finds the following equations of motion for the metric g pv 
and the antisymmetric tensor field B^ v 

\Bl p (B*B) 



B 



m 



B 



0. 



b, v - \b; u (b*b) 



G 



fJ,V 



kg, u + - ( 3H, Xp H/P - -g, u H 2 ) + 



mlB^-g^l + lB 2 



(2.6) 



(2.7) 



(2.8) 



1 + \ B2 - M B * B ) 2 

where V M is the covariant derivative. The equation (j2.7J) is also required from the con- 
sistency of (12. as you can easily check by applying V u to ()2.(ij) . Out of six equations 
in (jZHJ), one can find that three of them are dynamical equations and the other three are 
constraints. Therefore, the field B pu contains only three degrees of freedom as anticipated 
from the discussion of (jl.lj) . 



3 Time-dependent Homogeneous Solutions 

In this section we investigate evolutions of homogeneous magnetic components of the 
antisymmetric tensor field Bij with and without gravity. In the flat spacetime composed of 
a D-brane we find an oscillating solution. In the B-dominated D-brane universe, isotropy 
is recovered despite of the initially-assumed anisotropy, whereas it is not without D-brane. 



3.1 Flat spacetime: Homogeneous magnetic 5-matter 



We begin this subsection by turning off the gravitational coupling on the D-brane. We 
now consider homogeneous configuration B pu (t) and look for its time evolution in flat 
spacetime with g pu = 77^. Time component of the equation ()2.7|) is trivially satisfied and 
spatial components become 



E — B(E B) 



1 — E 2 + B 2 



E-B 



0. 



(3.1) 



where B i0 = (E) 1 and % = e 0ij k(B 
component of (|2.fij) reads 



) k . This equation (jH.lj) combined with the Oi- 
E = B(E B) . (3.2) 
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This allows a non-trivial solution with E = and B(£) 7^ 0. Then, the dynamical equation 
of the B field is summarized by 

d 2 B = -ml . (3.3) 

Integration of the equation (|3.3j) gives 

£ = -^ T B 2 + V cS (B), (3.4) 
2m B 



where £ is an integration constant, 14ff(B) = y/1 + B 2 , and overdot is differentiation over 
time variable t. Let us assume that the direction of the magnetic field is fixed, i.e., k 
defined by B(£) = B(t)k is independent of time. Then (13.4)1 is rewritten by an integral 
equation 

/ x f dB 

m B (t-t )=± ■ (3.5) 



JB ° ^2(£-VTTB^) 

The only pattern of nontrivial solutions is oscillating one with amplitude \/£ 2 — 1 for 
£ > 1. This oscillating solution with fixed amplitude is natural at classical level in flat 
spacetime when a homogeneous condensation of magnetic field is given by an initial con- 
dition. In expanding universes, the solution is expected to change to that with oscillation 
and damping due to growing of the spatial scale factor, which means the dilution of 
.B-matter. 



3.2 B-dominated universe 

We are interested in the cosmological implications of the large scale fluxes of the anti- 
symmetric tensor field, which might have existed in the early universe. Thus we look for 
cosmological solutions of (j2.6|) — f|2.8j) with appropriate initial conditions and examine the 
effects on the cosmological evolution. First we consider the case that the universe is dom- 
inated by the massive tensor field. Extending the above flat spacetime solution, we take 
into account the homogeneous configurations with an ansatz B^ = Bij(t) and B 0i = 0. For 
general non- vanishing Bij(t), the right-hand side of (J2.8)) says that the energy-momentum 
tensor of antisymmetric tensor field will have non-vanishing off-diagonal components 
Tij. To consider the simplest form of anisotropic cosmology, we take only single com- 
ponent of Bij to be nonzero; namely Bi 2 {t) = B{t) and B 2 3(t) = B 3 i(t) = 0. Even in 
this case, Tu(t) = T 22 (t) 7^ T 33 (t) in general and thereby the corresponding metric is not 
isotropic. Therefore, we take the metric to be of Bianchi type I 

ds 2 = -dt 2 + axitfidx 1 ) 2 + a 2 (t) 2 (dx 2 ) 2 + a 3 (t) 2 (dx 3 ) 2 . (3.6) 
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The energy-momentum tensor of the B field is written in the form 

(Tb)^ = (Tbulk)// + (Tbrane)// (3.7) 

= «7/ 2 A<J/ + diag [-p B , -p B , -p B , +Pb] + diag [-p b , -p b , -p b , -p b ] , 

where ps = B 2 / 'A/^lafal, p b = {m B /2K 2 ^)^Jl + B 2 / a\a\ and p b = (m 2 B /2Kl)/ \/l + B 2 jd\a^ 
Thus, the Einstein equations ()2.8|) for the tensor-field-dominated case are 



aia 2 a 2 a 3 a 3 ai 2 2 

1 1 = &\Pb + i^lPb + A, (3.8) 

Oi O2 02 o 3 3 oi 

5 2 S 3 S 2 S 3 2 2 ~ 1 a /o n\ 

1 1 = kIPb + KlPb + A, (3.9) 

Q-2 O3 O2 O3 

5 3 Si S3 Ol 2 1 2 ~ 1 A /QiM 

1 1 = k-IPb + K>tPb + A, (3.10) 

a 3 ai 03 ai 

Si S 2 Si S 2 2 1 2 1 A /o 1 1 \ 

1 1 = -K 4 p B + K 4 p b + A. (3.11) 



d\ CL 2 Oi O2 

The equation of motion for B(t) is 



' 2 fl 2 



01 a 2 o 3 y a/1 + fi 2 /af( 

Introducing = lnoj, we rewrite the Einstein equations f)3.8j) — f)3. 1 lj) as 

did 2 + d 2 d 3 + d 3 di = K 2 pB + n 2 p b + A, (3.13) 
ai + di(di + d 2 + 0:3) = «4P& + A, (3.14) 

d 2 + d 2 (di + d 2 + d 3 ) = K4P6 + A, (3.15) 

d 3 + d 3 (di + 0:2 + d 3 ) = 2K 2 p B + k 2 ^ + A. (3.16) 

Subtraction of ()3.15|) from ()3. 14|) gives 

di - d 2 + (di - d 2 )(di + d 2 + d 3 ) = 0. (3.17) 

A natural solution is di = a 2 which, with the help of scaling of x 1 and x 2 coordinates, 
leads to the isotropy in a; 1 x 2 -plane, ot\ = a 2 or equivalently ai = a 2 . 

3.2.1 Massless limit (m B = 0) 

To see the effect of the brane on the spacetime dynamics, let us first consider the limit 
of tub = which corresponds to either the absence of the brane or the limit of vanishing 
string coupling g s , namely the usual massless antisymmetric tensor field. In this limit, 
the equation ()3.12j) is easily integrated to yield a constant of motion 

— - — = L 3 (constant). (3.18) 

Oi0 2 
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With the vanishing potential, B manifests itself by non-vanishing time derivatives. In 
the dual variable, it corresponds to the homogeneous gradient along x 3 -direction. The 
spacetime evolution with the dilaton rolling in this case was studied in Ref. jHj- Here we 
have assumed that the dilaton is stabilized by some mechanism. Following Ref. jH] we 
introduce a new time coordinate A via the relation dX = L 3 dt j a\a 2 a 3 . Then the equations 
(I3.13|) - (j3.16|) can be written as 



a x a 2 + «2 a 3 + a 3 a i 
ol[ = a 2 

a" 



4 
0. 
1 



2 2 

-a-, a, 



1"2> 



2„2 



-a, a: 



1"2 ' 



(3.19) 
(3.20) 
(3.21) 



where the prime denotes the differentiation with respect to A. 
The solutions for ot\ and a 2 are trivial 

«i = CiA, a 2 = C 2 X, 



(3.22) 



where Ci i2 are constants and we omitted the integration constants corresponding simply 
to re-scaling of scale factors. The adequation (|3.21|) is also easily integrated to give 



3 2(Ci+C 2 )A 



a 3 



C 3 X. 



8(Ci + c 2 y 

The constraint equation ()3.19|) relates Ci j2 and C 3 by C 3 
relation between A and L 3 t is explicitly given by 

,A 

L 3 t — I dX ai(A)a 2 (A)a 3 (A) 

r-A 



(3.23) 

-C 1 C 2 /(C 1 + C 2 ). Then the 



dX exp 



(d + C 2 + C 3 )X + 



e 2(C 1 +C 2 )\ 



[8(d + C 2 )]^r 



(3.24) 



where x = e 2 ( Cl+C2 > A /8(Ci + C 2 ) 2 and P = (C x + C 2 - C 3 )/2(C 1 + C 2 ). The evolution of 
scale factors for large L 3 t is given by 



a 1)2 oc (log L 3 t) 



91,2 



a 3 oc L 3 t, 



(3.25) 



where gi j2 = Ci j2 /2(Ci + C 2 ). Therefore, with non-vanishing B 12 (t), only a 3 grows signif- 
icantly and the spatial anisotropy develops. This can be seen clearly by considering the 
ratio H 3 /Hi y2 where Hi = dj/aj. Taking C\ = C 2 and thus gi )2 = 1/4, we get 

#3 



Hi, 



41oeL 3 t 



(3.26) 
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which grows as time elapses. Finally, we remark that such anisotropy cannot be overcome 
by some other type of isotropic energy density in an expanding universe. Assume the 
isotropic universe (a, = a) driven by, e.g, radiation energy density pr. Then, one finds 
Pb oc 1/a 2 from (J3.7)l and (j3.18J) and thus Pb/pr oc a 2 , which implies that the late-time 
isotropic solution can be realized only in a contracting universe |6 a . 

3.2.2 The B oscillation 

Let us now take into consideration the effect of space-filling D-brane. To get sensible 
solution, we fine-tune the bulk cosmological constant term to cancel the brane tension, 
that is A = —m B /2, so that the effective four- dimensional cosmological constant vanishes. 
It is convenient to define a new variable b = B/a\a2- Then we can rewrite the full 
equations as follows 



ni 



+ 2di& 

(yi+& - 1 



b + 2a 1 b) + 



a\ + 2did3 
&\ + di (2di + 0:3) 
a 3 + a 3 (2di + d 3 ) 

b + (2ai + d 3 ) b + ^2«x + 2d 1 d 3 + 
With the help of (|3~27j) . can be written as 

b + (2di + d 3 ) b + 



mi 



2 \VT+¥ 



m 



B 



0. 



Aa\ + m 2 B Vl + b 2 - 1 + 



1 



0. 



(3.27) 
(3.28) 
(3.29) 

(3.30) 

VTTVJl- ■ (3 ' 31) 

which is valid only for the -B-dominated case, but more convenient for numerical analysis. 
Here we have put «4 = 1. 

First we examine the evolution of b(t) and scale factors qualitatively. From ([3.27)1 - 
(13.30)1 . we see that the natural time scale is m^ 1 and m^t is the dimensionless time 
variable. Suppose b starts to roll from an initial value bo, while the universe is isotropic 
in the sense that dio = 0:30- We assume initially a w = a 30 = 1 (ctio = am = 0) and 
B = so that b = B Q and b = — 2d 10 -Bo- While b is much larger than unity, the 
rapid expansion of cti due to the large potential proportional to m 2 B b drives b in feedback 
to drop very quickly to a small value of order one. Our numerical analysis in Figure 1 
shows that this happens within < 2 up to reasonably large value of 60 for which the 
numerical solution is working. The behavior of b(t) after this point is almost universal 
irrespective of the initial value 60 if h is much larger than unity. 

Once b becomes smaller than unity, the quadratic term of mass dominates over the 
expansion and b begins to oscillate about b = 0. Then the expansion of the universe 
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P3 



provides the slow decrease of the oscillation amplitude. The situation is the same as that 
of the coherently oscillating scalar field such as the axion or the moduli in the expanding 
universe. For small b, the energy-momentum tensor of the oscillating B field is given by 
T/ = diag[-p,pi,p 2 ,P3] where 

p = - (b + 2di&) 2 + l -m\ (Vl + V - l) « ~ (& 2 + m|& 2 ) , (3.32) 
p 1=P2 = -i(6 + 2d 1 6) 2 -i<(- 7 ^-l)^-i(6 2 -^), (3.33) 

J (b + 2d x 6) 2 - (v/fTfe 2 - l) « i (ft 2 - m|& 2 ) . (3.34) 

With the expansion of the universe neglected, the equation of motion for b is then ap- 
proximated by 

b + m 2 B b^0. (3.35) 

Since the oscillation is much faster than the expansion, we can use the time-averaged 
quantities over one period of oscillation for the evolution of spacetime. The equation 
(J3.35)) gives the relation (b 2 ) = (m 2 B b 2 ). Thus, the oscillating B field has the property 
PijP2,P3 ~ and behaves like homogeneous and isotropic matter. This justifies the name 
of B-matter. Therefore, after b begins to oscillate, the isotropy of the universe is recovered. 
To quantify how the universe recovers isotropy, let us define the parameter 

s = V2- ] 



2H X + H 3 

When Hi ~ H3, the evolution of the quantity s(t) is determined by 

1 (Api-pz- 9H 2 \ _ , lft-p 3 , ^ 2> 



where H = J2i H/3. At late time, the equation ()3.36|) is approximated as s ~ —(3/2)Hs 
and thus one finds s oc 1/t. Since H is also proportional to 1/t, the initial and final 
anisotropy Sij follow the relation 

H f 

s f = Si -f, (3.37) 

where Hij denotes the initial and final Hubble parameter, respectively. To get an idea of 
how fast the anisotropy disappears, let us consider Hf ~ 10~ 15 GeV which corresponds 
to the Hubble parameter at the electroweak symmetry breaking scale of temperature, 
T ~ 100 GeV. Taking the initial condition Sj ~ 1 around the beginning of the B oscillation 
Hi ~ tub, we find that the final anisotropy can be completely neglected for reasonable 
value of ttib- 
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For m B i> 1 with Hi ~ H 2 ^ H 3 and 5 1, the asymptotic solution of ()3.27|) - ()3.30|) 
can be explicitly found to yield 



#oc2/3t, bocl/t, p B , b otl/t 2 



(3.38) 



which shows the usual matter-like evolution. That is, the total energy density of the 
5-matter diminishes like p oc 1/a 3 even though the amplitude B(t) itself grows like 



Our qualitative results can be confirmed by solving the equations ()3.28j) . ()3.29j) . and 
(13.31)1 numerically. Figure 1 shows the numerical solutions for the initial value bo = 100. 
For the large value of bo, a± grows very fast while a 3 is frozen until b becomes smaller 
than unity. Then b begins to oscillate and the universe becomes isotropic again in that 
the expansion rates, di and a 3 , converge and finally the universe becomes B-matter 
dominated. 

3.3 Evolution of B(t) in the RW background 

When the universe is dominated by homogeneous and isotropic matter or radiation, while 
the B field contributes a minor fraction, the background geometry is described by the 
Robertson- Walker metric and the equation for b(t) is 



In this subsection, we study the evolution of b(t) during inflation and the radiation dom- 
inated era. 

3.3.1 Evolution during inflation 

During inflation, a/a ^ H (constant) and a/a ~ H 2 . We are assuming that the inflaton 
dominates over the antisymmetric tensor field, which means H 2 ^> m 2 B (\/l + b 2 — 1). 
Defining a dimensionless time variable as Ht, we write the equation (j3.39|) as 



Due to 2H 2 term, dominating over the oscillatory m 2 B term, b exponentially dies away. 
Neglecting m 2 B term, we have the solution b w boe~ 2Ht , corresponding to a frozen B 
field at some value B , while its energy density decreases exponentially during inflation. 
Therefore inflation dilutes away the pre-existing B field. 



B{t) oc a 



1/2 




(3.39) 




(3.40) 
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5 10 15 20 

m b t 

Figure 1: Numerical solutions for bo = 100. (a) The evolution of b(t). (b) Scale factors 
a>i(t) (solid line) and a 3 (t) (dotted line), (c) Expansion rates Hi(t) (solid line) and H 3 (t) 
(dotted line). 
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3.3.2 Evolution during radiation dominated era 

During the radiation dominated era, a/ a = l/2t and a/a = — l/4t 2 . Defining the dimen- 
sionless time variable as m B t, we rewrite the equation ()3.39|) as 

b + — 6+| , 1 r= + ^=L= 1 6 = 0. (3.41) 

As in the case of B domination, the amplitude of b shortly diminishes and becomes smaller 
than unity. Then the equation (|3.41j) is approximated as a linear equation in 6 to yield 
the usual damped-oscillation; 



bit) 2 



C\ cos I m B t + — J + C 2 sin (m B t + — 



7T (m B t) 3/4 

for mst ^> 1. In the regime, the energy density of the 5- matter is found to be 



(3.42) 



p oc j^j- 2 oc 1/a 3 . (3.43) 

That is, for mst ^> 1, the oscillation dominates over the expansion and the oscillating B 
field behaves like the usual matter during the radiation dominated era. As a consequence 
the radiation will be overtaken by the S-matter eventually and the cosmological evolution 
follows then the result of the subsection 3.2. 



4 Conclusion 

We have investigated the cosmology with the antisymmetric tensor field in a D-brane 
universe. A peculiar feature of this system is that the antisymmetric tensor field becomes 
massive due to its coupling to the massless U(l) gauge field on the brane and thus its 
cosmological consequences are drastically different from those of the massless antisymmet- 
ric tensor field. Analyzing the Bianchi type I cosmology describing the simplest form of 
an anisotropic universe, we found time-dependent homogeneous solutions exhibiting the 
matter- like (pressureless) behavior of the massive antisymmetric tensor field (5-matter). 
For instance, the 5-matter amplitude oscillating with the frequency of gets damped 
due to expansion. The energy density scales like p oc 1/a 3 and the averaged pressure 
becomes vanishingly small. As a consequence, the initial anisotropy dies away rapidly 
and the isotropic universe is recovered at late time. 
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